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Recent developments and results on the comparison of gluon to quark jets are discussed. A most important topic is the introduction 
of explicit energy scales of the jets. The scaling violation of the fragmentation function and the increase of the multiplicity with 
scale is shown to be directly proportional to the corresponding gluon or quark colour factor. The ratio of the hadron multiplicity in 
gluon to quark jets is understood to be smaller than the colour factor ratio due to differences in the fragmentation of the leading 
quark or gluon. Novel algorithms to reconstruct the colour portraits or the colour flow of an event are presented. 
Talk given at the XXIX ICHEP, Vancouver, July 1998. A shortened version will appear in the conference proceedings. 



1 Introduction 

QCD contains two types of colour charged fields, quarks 
and gluons. The relative coupling strength of the quark 
gluon vertex and of the triple gluon vertex is determined 
by the QCD colour factors Cf and Ca- The values of 
the colour factors directly follow from the SU(3) group 
structure of QCD. Reversely a precise measurement of 
the ratio of the couplings uniquely fixes the gauge group. 

A basic prediction thus is that there should be about 
Ca/Cf = 2.25 times mote gluon bremsstrahlung in gluon 
compared to quark_iets a. Scaling violation of the frag- 
mentation functionH'Q and the production of soft hadrons 
□ is correspondingly stronger in gluon jets. 

The comparison of gluon and quark jets thus allows 
for basic and intuitive QCD tests. Moreover this com- 
parison offers a direct handle on the fragmentation, the 
transition of colour charged fields into hadrons indepen- 
dent of fragmentation models. 

Due to colour confinement there are, however, no free 
quarks and gluons. The definition of gluon and quark 
jets therefore relies on the analogy to tree level graphs. 
Corrections of order a s /n have always to be expected. 
As QCD is a quantum theory also interference effects are 
anticipated due to the underlying event structure and 
due to the coherence of soft radiation. Both effects have 
indeed been observed. Finally effects due to the hadronic 
final state should be present. This is, however, part of 
fragmentation and what is to be understood. A final 
point of complication are the ambiguities introduced by 
the algorithms which are used to define the jets. 

2 Experimental Topics and Energy Scales 

The success of gluon and quark jet comparisons at LEP 
is due to two topics - the possibility of an unbiased tag of 
gluon jets and the better understanding of the (energy) 
scales underlying jet evolution. 

Gluon jets were originally identified in (2- or 3-fold) 
symmetric three-jet ( Y or Mercedes ) events where two 
jets could be identified as (heavy) quark jets using im- 
pact parameter tags or energy ordering (see e.g. u). The 
remaining jet is then a gluon jet. A similar light quark 



result can be obtained from a mixed quark/gluon sample 
from symmetric events by subtracting the gluon. Re- 
cently this technique -has also been extended to non- 
symmetric topologies Q. In this case one relies on the 
quark/gluon composition as predicted by the three-jet 
matrix element. The validity of this technique can, 
however, be tested by comparing results obtained with 
symmetric and non-symmetric events. This technique 
strongly improves the available statistic of gluon jets giv- 
ing access to a much wider range of energy scales and 
allows so to compare the dynamics of gluon and quark 
jets. 

A proper gluon to quark comparison has to be done 
with jets of a comparable scale. It turns out that this 
scale is not just the jet energy. The phase space for soft 
radiation is limited by coherence effects to cones given 
by the opening angles between the jets (angular order- 
ing). The relevant scale thus is in general a product of 
jet energy E times opening angle 9 i.e. a transverse mo- 
mentum. This is what has originally been used in the so 
called MLLA (modified leading log approximation) calcu- 
lations which were limited to small angles. If also larger 
opening angles are included it turns out that the so called 
hardness re — 2Esm.9/2 is a better choiceLra. This defini- 
tion corresponds to an equivalent CMS energy. It should 
however be kept in mind that in multi-jet events several 
scales are relevant, in principle. In so far the usage of re is 
an approximation which needs to be verified with data. 

3 Results 

3.1 Properties of Gluons at Fixed Scale 

So far measurements observed-gluon to quark multiplicity 
ratios in the range 1. to 1.5 EL This low ratio has been 
ascribed to missing energy conservation in the theoretical 
calculations Ej. To minimise this effect as welLas overlap 
effects between the jets it has been suggested til to study 
the most energetic gluons accessible at LEP, in events 
where the gluon recoils with respect t o th e qq -pair. The 
study of these rare events by OPAL t3 yields only 546 
gluons from the total data set taken at the Z resonance. 
The average gluon energy is «42GeV. The assignment 
of hadrons to the gluon is made inclusively: all hadrons 
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Figure 1: Gluon and quark fragmentation functions to charged hadrons as function of the scaled hadron energy ig (a) and the rapidity y 
(c) with respect to the sphericity axis as measured by OPAL. The lower plots (b,d) show the gluon to quark ratio. 



Table 1: Parameters of the gluon and quark multiplicity distribu- 
tions measured by OPAL. (iV c /,) is the mean charged multiplicity, 
D is the dispersion, 7 is the skewness and c is the kurtosis. 





Gluon 


Quark 


(Nch) 
D 

7 
c 


14.32 ± 0.23 ± 0.40 
4.37 ± 0.19 ± 0.26 
0.38 ± 0.13 ± 0.18 
0.18 ± 0.34 ± 0.30 


10.10 ± 0.01 ± 0.18 
4.30 ± 0.01 ± 0.10 
0.82 ± 0.01 ± 0.04 
0.98 ± 0.03 ± 0.11 



in the gluon hemisphere of the event are assigned to the 
gluon. 

This selection proves to be rohiist against the us- 
age of different cluster algorithms El A comparison of 
these inclusive gluon-jets with identified uds quark jets 
yields multiplicity distributions with parameters given in 
Tab. [|. Skewness, kurtosis and multiplicity differ signifi- 
cantly between the gluon and quark distributions. If the 
slightly higher energy of the quark jets is considered by a 
Monte Carlo correction, even for these high energy glu- 
ons the multiplicity ratio still attains, values much smaller 
than the simple QCD expectationl!3: 



N 



gluon 
ch 



N 



quark 



= 1.509 ± 0.022 ±0.046 



(1) 



A comparison of the gluon and quark factorial and 
cumulant moments has also been performed by OPAL. 
Cumulants K q of order q measure the genuine g-particle 



correlation. The predictions of these quantities w< 
only found to agree "qualitatively" with the data 
Partly this seems understandable as the leading order 
(LO) , next to LO (NLO) and NNLO predictions dif- 
fer strongly__A recent numerical calculation 113 continuing 
the work ofliS, which now also considers energy conserva- 
tion, the correct low energy behaviour as well as OPAL's 
inclusive gluon jet definition, describes the data well. 

Also charged hadron spectra have been studied for 
the inclusive gluonstLp, see Fig. [l]. Consistent with previ- 
ous studies (see e.g. □) a roughly exponential falloff is ob- 
served for the distribution of the scaled charged hadron 
energy, xe, which is more pronounced for gluons than 
for quarks. At large Xe the gluon distribution is sup- 
pressed with respect to the quark result by more than 
one order of magnitude, i.e. the production of fast, so- 
called leading particles is strongly suppressed in gluon 
jets. At the smallest xe, however, the gluon to quark 
ratio reaches values of up to 1.8. The distribution of the 
rapidity with respect to the sphericity axis of the event 
shows a corresponding behaviour. This choice of the ob- 
servable "zooms" the range of low hadron energies. At 
the smallest rapidities (y <j-3) a high gluon to quark ratio 
is observed as expected inO: 

dN sluon /dy 



dN% ar Vdy 



= 1.815 ±0.038 ±0.062 



(2) 



Although not mentioned mill it is important to note that 
in the rapidity range 3 < y < 6 the gluon distribution 
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Figure 2: Gluon and quark splitting kernels (left) and their ratio (right) as measured by DELPHI. Theoretical expectations are shown as 
lines. 



is below the quark distribution. The area between the 
gluon and the quark distribution amounts to almost one 
unit of charged multiplicity. This suppression further 
extends down to rapidities of y ~ 2 as is observed from 
Fig. |l]d) and necessarily needs to be taken into account 
when comparing the QCD expectation for the gluon to 
quark multiplicity ratio to data, e.g. for Eqn. ^. 

In order to verify the QCD prediction of a stronger 
splitting probability for gluons and also to clarify the 
origin of the small observed multiplicity ratio DELPHI 
directly measured the so-called splitting kernels. These 
are the most important piece of the evolution equations 
for the generating functions or functionals which describe 
the jet- and hadron multiplicities (see e.g. El). The im- 
portance of these kernels, commonly known is the special 
(collinear) case of the Altarelli-Parisii kernels, can be un.- 
derstood from the following simplified argumentation □. 
Consider a given number of partons Ni (y) which have not 
split (or decayed) at a given jet resolution y. The num- 
ber of partons which will split in an interval Ay then is 
expected to be: 



AN^y) = -F(y) 
1 



f{v) = 



N 1 (y)-Ay 
dNi(y) 



(3) 



Ni(y) dy 

This equation is analogous to the radiative decay equa- 
tion, however, with a y dependent decay "parameter" 
F(y) describing the dynamics of the process. This is the 
splitting kernel. Eqn. || can be expressed by the usual 
jet rates R n — N n /N tot and then yields: 



F(y) = D 1 (y) 



rflog Rxjy) 
dy 



(4) 



In analogy to the differential jet rates known form the a s 
analyses, the quantity D\(y) is called modified differen- 
tial 1-jet rate. It is a direct experimental measure of the 
splitting kernel. 
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Figure 3: Average subjet multiplicity minus 1 for gluon and quark 
jets (upper plot) and gluon to quark ratio (lower plot) as function 
of the resolution y. 
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Figure 4: Left side: Scale dependence of the quark fragmentation function measured in three-jet events as function of the hardness scale 
K compared to measurements in e + e - annihilation. Lines are power fits to the three-jet data. Right side: Scale dependence of the gluon 
fragmentation function measured in 3 jet events as function of the hardness scale k. Data points belong to symmetric and non-symmetric 
topologies. Lines result from a simultaneous DGLAP fit of the gluon and quark fragmentation function. 



The result of the DELPHI measurement for gluons 
and quarks as function of y as defined with the Durham 
algorithm B as well as the corresponding ratio is shown 
in Fig. ||. At large y the data about follow the perturba- 
tive expectation. A deviation at very large y (best seen 
in the gluon to quark ratio) is .-due to the structure of 
the underlying three-jet events a (see also behaviour of 
the subjet multiplicity- 1 inH). The gluon and the quark 
cannot be treated as independent objects in this limit as 
the large y splitting resolves the underlying event. At 
smaller y both the quark and the gluon data fall below 
the perturbative expectation. This is to be expected due 
to non-perturbative hadronization effects. It is most im- 
portant, that this influence sets in already at larger y (i.e. 
"earlier" ) for gluons compared to quarks. This causes the 
initially large gluon to quark ratio to fall off rapidly with 
falling y. 

In the framework of the HERWIG cluster fragmenta- 
tion modeltj this effect can easily be understood. Quarks 
are valence particles of hadrons whereas gluons are not. 
Gluons have to be split into a qq -pair first. So hadroni- 
sation effects should be visible already at larger y in case 
of gluons. The same effect is responsible for the suppres- 
sion of the fragmentation function at large xe or large 
rapidity y as discussed above. 

The behaviour of the so called sub-jet multiplicity 
< rij e t >= Yli i ' Ri ( see Fig- H 8 -)) illustrates that this is 
indeed the reason for the gluon to quark multiplicity ratio 
to be smaller than expected. Initially at logy w —2 the 
subjet multiplicity- 1 ratio [] (see Fig]|b) reaches about 
the expected value ~2. Then in coincidence with the 

a l is subtracted to consider the initially present parton. 



deviation of Df (y) from the perturbative expectation this 
ratio falls off and rapidly reaches values typical for the 
hadronic multiplicity ratio. Note from Fig. ||a) that this 
happens when the typical subjet multiplicities are around 
3, thus at much larger y than most hadrons are formed. 

3.2 Scale Dependence of Quark and Gluon Fragmenta- 
tion 

DELPHI in several publicationsii : 0i attempted to com- 
pare the scale dependence thus the dynamical evolution 
of gluons and quarks. Using the full data set containing 
about 85000 identified gluon jets a detailed comparison 
of the gluon and quark fragmentation function D{xe) as 
function of k is performed a. The measurement is shown 
in Fig. ^. The quark fragmentation function as obtained 
from quark jets in three-jet events (left side) was fitted by 
a power ansatz superimposed to the data (open points). 
Also shown (full points) is the fragmentation function 
as measured in e + e~ annihilation from PETRA energies 
E2l up to the highest LEP energies. Both measurements 
agree very well, enforcing the interpretation of k as a 
valid scale for jet evolution in three-jet events. It should, 
however, be mentioned that the transverse momentum of 
a jet with respect to the event axis works similarly well. 
The corresponding result for gluons is shown on the right 
side of Fig. [|. Results for symmetric and for all topolo- 
gies agree well, strengthening again the validity of the 
analysis and the scale. The fits shown with the scale de- 
pendence of the gluon fragmentation function result from 
a simultaneous 1 st order DGLAP evolution together with 
the quark fragmentation function. The strong coupling, 
the colour factor ratio Ca/Cf, as well as the parameters 
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Figure 5: Slopes of the scale dependence of the quark and gluon 
fragmentation function as function of xe- The shaded area indi- 
cates the stat. uncertainty of Ca from the DGLAP fit. Data points 
result from power law fits of the individual xe intervals. 



of the analytic ansatz of the fragmentation functions were 
treated as free parameters. The fit yields the following 
preliminary result for the colour factor ratio: 



Ca 
C f 



2.44 ± 0.21 (stat. 



(5) 



Fig. H compares the logarithmic slope of the gluon 
and quark fragmentation function as function of xe- The 
function is the result of the DGLAP fit of the data. The 
shaded area indicates the uncertainty on Ca for fixed 
Cf = 4/3. The data points were obtained from power 
law fits to each gluon and quark xe interval individ- 
ually. As expected from the structure of the DGLAP 
equation scaling violations are much stronger for gluons 
than for quarks [| Due to the stronger scaling violations 
in gluon jets pp sscattering which is dominated by gluon- 
gluon scattering should be well suited to measure a s form 
the scaling violation of the charged hadron fragmentation 
function. 

It is evident from the gluon to quark comparison at 
low xe - indicating a ratio of about 2 (see Fig. ||) - that 
the scale dependence of the hadron multiplicity opens a 
new way in testing the prediction: 



N' 



gluon 



ch 



ch 



Ca 
C f 



(G) 



insensitive to the non-perturbative effects discussed 
above. Note that the prediction Eqn. ||is expected to be 
valid for any large scale if non-perturbative effects can 

b Note that the Altarelli-Parisii splitting functions include the in- 
dividual colour factors multiplicatively. 



be neglected. Consequently also the ratio of the deriva- 
tive of the multiplicities with respect to the relevant scale 
(here taken to be k) is equal to the colour factor ratio: 



8 < Nf h uon (K) > /dn 



Ca 
C f 



(7) 



This prediction is equivalent to Eqn. || but, as the mul- 
tiplicity evolves slowly with scale, should hold already at 
much smaller scale. The comparison of the gluon and 
quark splitting kernels as well as the behaviour of the 
fragmentation function suggest that most of the non- 
perturbative difference between gluons and quarks is due 
to the fragmentation of the leading quark or gluon, thus 
is to be expected to happen at very small scale. A cor- 
responding calculation in the framework of the colour 
dipole modeled confirms this expectation. 

Fig. H compares the dependence of the multiplicity 
assigned to gluon and quark jets as function of log k. It is 
immediately evident from this plot that the slope of the 
gluons is roughly twice as big as for quarks thus already 
impressingly confirming the QCD expectation Eqn. 0. 
The data was fit with the predictions for—the evolution 
of the multiplicity with scale as given inc3E3, however, 
adding constant terms to account for non-perturbative 
differences between gluons and quarks, thus: 



< N > q ch («) 



<N>l h (K)=<N>H 



ch 

perturb. 



(«) 



N° (8) 



g q(n) is-a-NNLO prediction for the multiplicity ratio as 
3. The ansatz represents the data well. The 



given m 
resulting colour factor ratio is: 



Ca 
C f 



2.U±0.10{stat.) 



(9) 



No systematic error is given with this measurement as 
a superior method is discussed also in □ which deter- 
mines the colour factor ratio from the scale dependence 
of the multiplicity of three-jet events. It is insensi- 
tive to ambiguities induced by the choice of the scale 
and also accounts for coherent radiation from a qqg - 
ensemble (see also E3) . The non-perturbative term in- 
troduced in this study which accounts for differences 
in the fragmentation of the primary quark or gluon is 
A 9 - A 9 = 1.97 ± 0.03 ± 0.33. 

It is instructive to estimate a lower limit for Nq — Ni. 
from the gluon and quark fragmentation functions □ 
in the momentum range where the gluon is below the 
quark fragmentation function (compare also Fig.^a,c)). 
This limit can be calculated by subtracting the observed 
gluon multiplicity from the (wrongly) expected one: 
N ex P ected _ 2 N q . This yields A 9 - A S > 1.22 ± 0.08(> 
1.15 ±0.2) from Y and Mercedes events, respectively. As 
the suppression persists to lower xe or y (see Fig. np,d)) 
it probably amounts for the full non-perturbative term 
N«-N°. 

The left side of Fig. g presents the gluon to quark 
multiplicity ratio, the ratio of the fitted functions Eqn. 




Figure 6: Left: Scale (k) dependence of the jet multiplicity in gluon and quark jets as measured by DELPHI. Lines are the result of a fit of 
Eqn. [?] Right: Ratio of the gluon to quark multiplicity as function of the scale. The lower line is the ratio of the fits on the right. The upper 
line is the ratio of the slopes of these fits. 



U as well as the corresponding slope ratio. Also included 
are measurements _jat very small scale (4-7 GeV) per- 
formed by CLEO c3 in T(1S) — > gg7 decays and the 
OPAL result for inclusive gluons. For this point k has 
been estimated from_the average gluon energy and the 
angle cuts given in E3. As this result compares gluon 
to pure light quark (uds) jets it may be expected to lie 
slightly higher than the extrapolation of the DELPHI 
result. The fit describes all measurements well, confirm- 
ing that the non-perturbative effects can be absorbed 
effectively into constant terms allocated to the leading 
particles in a quark or gluon jet. 

3.3 Reconstruction of Colour Connections 

The stronger radiation off a gluon may also be inter- 
preted as due to two colour connections compared to one 
in the case of quarks. One may then try to reconstruct 
the colour connections in an event to distinguish e.g. 4 
quark final states (like WW — > qqqq events, 2 connec- 
tions) from QCD background (qqgg, 3 connections). One 
may also try to reduce the combinatorial background in 
combining only the colour connected qq -pairs in WW 
events. Besides several theoretical proposals related to 
this .subject 123 now first experimental results are avail- 
ablec§E3. These are based on the idea of reconstructing 
the colour portrait of a processE3. To introduce the algo- 
rithm recall first the structure of the soft hadron (labeled 
h) production cmss^section in a hard n parton (labeled 
initial stateEa'B 



da 



soft 



da ha 



■ dn 



dEh 
E h 



n 

i,j hard 



c« • if. 



= sm%/2 
13 sin 2 6 lh /2 ■ sin 2 6 jh /2 



(10) 



Important pieces in this equation are the colour factors 
Cij which contain the information about the colour con- 
nections (Cij — if no connection) and the so called 
antenna pattern Wij which determine the angular struc- 
ture of the radiation. The poles O^h ~ > appear in the 
direction of the hard partons and lead to the formation of 
the jets. The nominator contains the angular dependence 
known from the hardness scale. The idea is to reconstruct 
the colour connections by adequately weighting the ob- 
served soft hadron production. The algorithm proceeds 
as follows: 

• Reconstruct jet directions from leading particles. 

• Calculate the probability for un-associated (low en- 
ergy, large angle) particles to stem from a cluster 
(hard parton) j: 



K hj 



C 3 



2E£ sin" hj /2 



(11) 



Renormalise so that Whj = 1 • 



• Associate hadron h to the cluster pair which max- 
imises Whu = u>hi + Whk ■ Calculate the colour con- 
nection coefficient: 



W kl =C Y, 9(E h ) 



■ Whkl 



(12) 



h assoc. 



g(Ei) is an energy dependent weight, C is chosen 
so to properly normalise the Wti] Wm = 2. 




Figure 7: Fish-eye plot of the hadron distribution in a typical 
Mercedes event (Monte Carlo). 
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Figure 8: Weight distribution for Mercedes events (DELPHI data). 



To test the technique it has been applied to Mercedes 
events (Monte Carlo and data) where the (heavy) quark 
jets wcre_additionally identified using an impact param- 
eter tagc3. Fig.fj] shows a typical event. It is clearly 
observed that the soft hadrons (boxes) accompany the 
quark-gluon colour connections indicated by the lines. 
The corresponding colour connection coefficients clearly 
distinguish the qg(qg) from the qq case. Fig. ||| (left) 
shows the colour connection coefficients which were ob- 
tained from DELPHI events of the Mercedes topology. 
The additional impact parameter identification allows to 
measure directly the purity and efficiency of the algo- 
rithm. The data show that purities of 90% at 80% effi- 
ciency can be obtained. 

A modification of the algorithm has also been ap- 
plied to the analysis of WW — > qqq'q' events at 189 
GeV centre-jQf-mass energy measured by DELPHI. As 
proposed inE3 all (hypothetical) strings connecting the 



4 parton final state were boosted to their rest systems. 
Then the transverse momentum of the hadrons with re- 
spect to the string was calculated and similar to Eqn. 
[Tl| weights were defined using the transverse momentum 
of the hadrons iv^j — Ch/Pt- The string configuration 
which minimises the sum of the transverse momenta was 
selected. Fig. |9| shows the mass spectrum of all possible 
jet-jet pairs (upper plot) and taking only the combina- 
tions selected by the algorithm outlined above. A clear 
reduction of the combinatorial background due to wrong 
qq' pairing is indicated. Contrary to kinematic cuts the 
method does not bias the W-mass. So far no attempt 
has been made to reject QCD background. Thus the 
appearend reduction of this contribution is due to the 
reduced combinatoric only. Currently the method can 
only be used as an additional tool as the efficiency is only 
about 50%. Improvements are however to be expected 
in the near future. 
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Figure 9: Jet-jet mass distribution for four-jet events measured 
by DELPHI at 189 GeV. The upper plot contains all combinations, 
the lower only those selected with the algorithm described in the 
text. Histograms indicate the Monte Carlo predictions for WW 
and QCD events. 



4 Summary 

The introduction of energy scales in the comparison of 
gluon and quark jets opened the possibility for studies of 
the scale dependence of gluon and quark jet properties. 
This in turn lead to the understanding of the (smaller 
than originally expected) gluon to quark multiplicity ra- 
tio due to non-perturbative effects in the fragmentation 
of the leading gluon or quark and to the possibility to 
really measure the colour factor ratio from a direct gluon 
to quark jet comparison. These observations will soon 
be further exploited in order to identify quark and gluon 
jets or colour connections in multi-jet environments. 
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